Electronic Structure of the Complex Hydride NaAlH 4 
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Density functional calculations of the electronic structure of the complex hydride NaAlEU and 
the reference systems NaH and AIH3 are reported. We find a substantially ionic electronic structure 
for NaAltU, which emphasizes the importance of solid state effects in this material. The relaxed hy- 
drogen positions in NaAlHU are in good agreement with recent experiment. The electronic structure 
of AIH3 is also ionic. Implications for the binding of complex hydrides are discussed. 
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The complex hydrides, A x (MH4) y , with A=Li, Na, structure of solid NaAlH 4 is better described as mixed 
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K, Mg, Ca, Sr or a mixture of these and M=B or Al 
seem promising for H storage since they contain very high 
weight percent H, and much of the H content can be 
evolved at moderate temperatures. 1-5 However, hydro- 
gen desorption from these materials is associated with 
their melts, and they were long viewed intrinsically non- 
cyclable. However, in 1997 Bogdanovic and Schwickardi 
reported that with certain metal additions, particularly 
Ti, NaAlH 4 can be cycled. 1 This discovery opened the 
door for consideration of these complex hydrides as cy- 
clable hydrogen storage materials. However, this result 
has not been reproduced for other related complex hy- 
drides and understanding of the basic materials prop- 
erties that govern the uptake of H is still not complete. 
What is known is that hydrogen desorption from NaAlH 4 
takes place via a two step process, forming first NaaAlHg, 
Al and hydrogen, and then NaH and Al metal. Desorp- 
tion is associated with the temperature where melting 
occurs. While the Ti addition is often referred to as a 
catalyst or dopant, the actual role of Ti in enabling the 
cyclability of NaAlH 4 has yet to be established. Pos- 
sibilities include catalysis as mentioned, modification of 
the thermodynamics of the decomposition, i.e. the bal- 
ance between solid state Ti containing NaAlH 4 and the 
decomposed NaH - Al - Ti mixture, modification of the 
microstructure of post decomposition NaH + Al mix- 
ture, e.g. by keeping the precipitated Al and the NaH 
from segregating over large distances, e.g. by enhanc- 
ing the solubility of Al in NaH or by the formation of 
Ti-Al alloys, like T1AI3, instead of precipitated Al, and 
others. 2 ' 3 ' 7 " 12 

Developing understanding of the various phases in- 
volved will likely be important for sorting out the physics 
of cyclable hydrogen storage in this material. Here we 
start with the simplest of these phases, NaAlH 4 , by using 
density functional calculations of the electronic structure, 
in comparison with results for NaH and AIH3. Conven- 
tionally, the bonding in NaAlH 4 is viewed as that of a salt 
made of Na + cations and AIH4 anions, with the internal 
bonding of the AIH4 units being primarily covalent, con- 
sistent with the tetrahedral coordination of Al and what 
is expected in the liquid. We show that the electronic 



ionic, i.e., Na+AP+Hj . Solid NaH is found to be ionic 
as expected. AIH3 is also ionic but with a smaller band 
gap. The ionic nature of solid NaAlH 4 is understood as 
a result of long range Coulomb interactions, implying a 
greater sensitivity of the electronic structure and there- 
fore bonding of H in NaAlH 4 to substitutions and defects, 
than would otherwise be the case. Further this provides 
an explanation for the association between the hydrogen 
desorption and melting. 

The present calculations were done within the local 
density approximation (LDA) to density functional the- 
ory, using the general potential linearized augmented 
planewave method with local orbital extensions, as im- 
plemented in the WIEN2K code. 13-15 For consistency, 
the same LAPW sphere radii of 1.7 eto and 1.1 ao were 
used for the metal and hydrogen atoms, respectively, in 
all three compounds. Well converged basis sets consist- 
ing of an LAPW cutoff, k max =5.91 Oq 1 plus local orbitals 
were used (the effective dimensionless values of the ba- 
sis cutoff were Rk max — 6.50 for H and Rk max = 10.05 
for the metal atoms. The Brillouin zone samplings were 
done using the special k-points method, with 21, 38 and 
84 points in the irreducible wedge for NaAlH 4 , AIH3 and 
NaH, respectively. This was found to be well converged 
for these insulating materials. 

NaAlH 4 occurs in a tetragonal structure (space group 
Mi/a) with lattice parameters a=5.021 A, c=11.346 A, 
and its own structure type, which has two formula units 
per primitive cell. 16 Recent neutron measurements for 
NaAlD 4 17 confirm this structure, with slightly lower lat- 
tice parameters, but rather different H positions. NaH 
occurs in the NaCl structure, with lattice parameter 
a=4.88 A. 18 AIH3 occurs in a rhombohcdral structure 
(space group R3c or possibly i?3), with hexagonal lattice 
parameters, et=4.451 A, c=11.766 A, and two formula 
units per primitive rhombohedral cell (from neutron and 
X-ray diffraction). 19-21 In general, the structures of hy- 
drides, particularly the H positions, are difficult to de- 
termine because the H has a very small X-ray scatter- 
ing factor and samples can differ in H stoichiometry and 
ordering. Here we fully relaxed the internal structures 
using LDA total energies and forces, keeping the pre- 
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sumably reliable lattice parameters fixed at the reported 
experimental values. For NaH there was no relaxation to 
do because both atoms are on high symmetry sites. For 
NaAlH4 we obtain positions differing from the positions 
reported in Rcf. 16, but in close agreement with very re- 
cent neutron results. 17 The relaxed H positions (Wycoff 
notation, site 16/) are x=0.2364, y=0.3906 and z=0.5451 
as compared to x= 0.2371, y=0.3867 and z=0.5454 from 
neutron scattering. This structure is illustrated in Fig. 
1, which clearly shows the A1H 4 building blocks. The Al- 
ii bond length in our structure is 1.652 A. This is only 
a little bigger than the sum of the covalcnt radii of Al 
and H (0.37A+ 1.18 A= 1.55A) and would seem to be 
a reasonable number for covalently bonded A1H 4 units. 
For AIH3 we relaxed in the lower symmetry R3 space- 
group considered by Zogal et al. Our positions differ 
somewhat from the X-ray structure of Ref. 21, but are 
consistent with the assignment of R3 in the absence of 
H disorder. The Al - H nearest neighbor distance in our 
structure is 1.731 A, which is longer than in NaAlH4 and 
is also longer than the sum of the Al and H covalcnt radii. 
It should be mentioned that this Al - H bond length is 
close to the value of 1.715 A, from the structure of Turlcy 
and Rinn. 19,22 We also considered R3, but obtained no 
further relaxation (see below). 

The calculated LDA band structure and corresponding 
electronic density of states (DOS) for NaAlH 4 is shown 
in Figs. 2 and 3, respectively. The total DOS is similar 
to that recently reported by Vajeeston and co-workers, 23 
and is also similar to that reported for the related com- 
pound LiAlELi. 24 The band structure has a large ^4 eV 
band gap, separating H derived valence bands from metal 
derived conduction bands. We emphasize that despite 
the seeming AIH4 units in the structure, and the ex- 
pected covalency of such chemical units, the calculated 
electronic structure is very strongly ionic. In particular, 
it can be seen from Fig. 3 that the valence bands are 
strongly dominated by H, while the conduction bands 
have very much less H character. The valence bands 
consist of two crystal field split manifolds, each ~3 eV in 
width. The calculated DOS of NaH and A1H 3 are shown 
in Figs. 4 and 5, respectively, along with the projections 
onto the H LAPW spheres. Since, the 1.1 a spheres are 
not large enough to fully contain the Is states of ions, 
the H projection underestimates that H contribution to 
the electronic structure. However, the ratio between the 
projections on H from different energy regions, gives a 
good indication of the ratio of the H contributions to the 
electronic structure in those energy regions. The band 
structure of NaH is also strongly ionic, with a band gap 
slightly smaller than that of NaAlH 4 , while AIH3 is also 
ionic, but has a smaller ^2 eV band gap. The valence 
band width of AIH3 is ^9 eV. It should be noted that 
Goncharenko and co-workers 20 , had already conjectured 
that AIH3 is ionically bonded based on its crystal struc- 
ture. 

The reason for the ionic electronic structure of NaAlH 4 
can be understood as due to the long range Coulomb in- 



teraction in solids. This Ewald contribution to the energy 
favors ionic electronic structures, and is well known to 
stabilize 2 ~ in metal oxides, for example, even though 
dimers and small molecules with the same metal - O 
neighbors may be covalent. Here H~ is stabilized in this 
way. In metal oxides, especially when there is some co- 
valency between the O and nominally unoccupied metal 
orbitals, the O 2- ions are highly polarizable, as may be 
expected from the fact that 2 ~ outside the Coulomb 
field of the solid is not a stable ion. Following the argu- 
ments of Cohen, 25 which related ferroelectricity in oxide 
perovskites to ionic electronic structures with weak cova- 
lency of this type, and considering the electronic struc- 
ture of AIH3, it seemed worthwhile to check if ferro- 
electricity is present. Accordingly, we made small dis- 
placements of the atoms away from the relaxed posi- 
tions within the reduced symmetry non-centrosymmctric 
spacegroup i?3 and calculated the restoring forces. How- 
ever, no ferroelectric instability was found within this 
symmetry. 

Returning to the bonding of NaAlH4, we note some 
expected consequences of the ionic electronic structure. 
First of all, since the bonding is stabilized by long range 
interactions, rather than primarily short range Al-H co- 
valent bonds, it should be more sensitive to stoichiome- 
try, defects, lattice parameter changes and off-site sub- 
stitutions, than in a salt made of strongly covalcnt A1H 4 
units. This implies tunability of the hydrogen binding, 
e.g. by alloying, which in turn would offer tunability of 
the thermodynamic balance between the solid and the 
dehydrided NaH + Al mixture. Secondly, it provides a 
natural explanation of why the hydrogen desorption is 
strongly connected with melting. Presumably, melting 
involves disruption of the H lattice and with it a loss 
of the long range Coulomb stabilization of the H~ ions. 
The resulting loss of binding at melting then would re- 
sult in H release from the material. Finally, we note that 
while our calculations are specific to NaAlH 4 , the DOS 
of LiAlH 4 is qualitatively similar, 24 , suggesting that sim- 
ilar physics may be operative there and in other related 
complex hydrides. 
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FIG. 1. Crystal structure of tetragonal NaAlH4 with the 
relaxed atomic positions. The small spheres are H, the large 
dark spheres are Na and the large light spheres are Al. 
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FIG. 2. LDA band structure of NaAlH4 with the relaxed 
crystal structure. The band gap is between a H derived va- 
lence band and metal derived conduction bands. 
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FIG. 3. Electronic density of states and projection onto the 
H LAPW spheres for NaAlHU, on a per formula unit basis. 
Note the ionic nature shown by the very different hydrogen 
contributions to the valence and conduction bands. 

2.0 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 




E(eV) 



FIG. 4. Electronic density of states and projection onto the 
H LAPW sphere for NaH. 
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FIG. 5. Electronic density of states and projection onto the 
H LAPW spheres for AIH3, on a per formula unit basis, using 
the relaxed crystal structure. Note the smaller band gap but 
still ionic nature. 
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